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We present a finite-size scaling analysis of the entanglement in a two-dimensional arrays of quan- 
tum dots modeled by the Hubbard Hamiltonian on a triangular lattice. Using multistage block 
renormalization group approach, we have found that there is an abrupt jump of the entanglement 
when a first-order quantum phase transition occurs. At the critical point, the entanglement is 
constant, independent of the block size. 
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O ■ I. INTRODUCTION 

o ; 

After the pioneering work of Benioff 0, Feynmanj^], and Deutsch Q in field of quantum information and computing, 
the development in this field has been explosive in the past few years 0, 0] . A lot of exciting progresses have been 
made, such as the experimental realization of Shor's quantum factoring algorithm with NMR scheme for quantum 
computation , teleportation Q , cryptography j| and dense coding Q . Central to all these remarkable achievements 
is the concept of the quantum entanglements O, IL? ! K$L Il-ij . Just like energy, the entanglement has now been 
regarded as a controllable and a fungible physical resource [l5j . Almost all the efficient protocols in quantum computing 
applications are based upon its generation. Experimentally, entanglements have already been produced among up to 
four photons |l6l Il7| and even between two macroscopic states such as two superconducting qubits, each of which 
contains as large as 10 9 electrons |l8|. But theoretically, because an ensemble's Hilbert space grows exponentially with 
the number of its component particles, we are still far from fully understanding the contents of the entanglements. 
Only for the simplest state with two distinguishable particles can we have a complete description of the entanglement 
measure. For states of more than two particles, especially for mixed states, the current knowledge about their 
f"^ ' entanglement is very limited and all the related complexities have just begun to be explored. 
■ 

o ■ 

II. ENTANGLEMENT MEASURE FOR FERMIONS 

a: 

, 1 * . For spin-only entanglement of localized distinguishable particles, the most suitable and famous measure of the 
entanglement is the Wootters' measure Recently, Schlieman[2(],El examined the influence of quantum statistics 
upon the definition of the entanglement. He discussed a two-fermion system with the Slater decomposition instead of 
Schmit decomposition for the entanglement measure. If we take each of the indistinguishable fcrmions to be in the 
single-particle Hilbert space C N with f m , fm( m = 1> denotes the fermionic annihilation and creation operator 

of single-particle states and |A) represents the vacuum state. Then a pure two-electron state can be written as 
X ^rnnfmfn where ui mn = —Lu nm . In analogy to the Schmidt decomposition, it can be proved that every |\&) 
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can be represented in a appropriately chosen basis in C N in a form of Slater decomposition [2C 

K 
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where , « | A) , ,~ | A) , i — 1 , • • • , K, form an orthonormal basis in C N . The number of the nonvanishing coefficients 
Zi is called the Slater rank, which is then used for the entanglement measure. With similar technique, the case of 
two-boson system is studied by Li and Paskauskas 

Gittings |24| put forward three desirable properties of any entanglement measure: (a) Invariance under local unitary 
transformations; (b) Non invariance under non-local unitary transformations; (c) Correct behavior as distingusisha- 
bility of the subsystems is lost. These requirements make the distinction between one-particle unitary transformation 
and one-site unitary transformations become relevant. As claimed in Ref. |24j. a natural way achieving this distinc- 
tion is to use a basis based upon sites rather on particles. Through Gittings' investigation, it is shown that all the 
above-discussed entanglement measures fail the tests of the three criteria. Only the Zanardi's measure [25| survives, 
which is given in Fock space as the von Neuman entropy, namely, 

Ej = -Tr Pj hxpj,pj = Trj \ip){ip\, (2) 



2 



where Ttj denotes the trace over all but the jth site and ip is the antisymmetric wave function of the studied 
system. Hence Ej actually describes the entanglement of the jth site with the remaining sites. This measure is well 
operational numerically and will be used in our lattice systems, which uses the basis based upon sites from the start. 
A generalization of this one-site entanglement is to define an entanglement between one L-site block with the rest of 
the systems [29j . 

E L = -Tr{ PL \og 2 p L ). (3) 
In this paper, this idea will be pursued together with the real-space renormalization group (RG) method. 

III. ENTANGLEMENT AND QUANTUM PHASE TRANSITION 

Quantum phase transition (QPT) is a qualitative ch ang e of some physical properties in a quantum many-body 
system as some parameter in the Hamiltonian is varied |26j . At the critical point when QPT happens, a long-range 
correlation can develop in the system. QPT is caused by quantum fluctuations at the absolute zero of temperature 
and is a pure quantum effect. Similarly entanglement is also a unique quantum property with non-local states of two 
or more quantum systems superposed with each other. Hence it is inviting to study the relationship between them. 
Recently it has been argued that the property responsible for the long-range correlation in QPT is entanglement [27| • 
So far, there have been some efforts in this direction, such as the analysis of the XY model about the single-spin 
entropies and two-spin quantum correlations j2^| , the entanglement between a block of L contiguous sites and the rest 
of the chain [2^| and also the scaling of entanglement near QPT 30] . But because there is still no analytical proof to 
validate the above speculation, the role played by the entanglement in quantum critical phenomena remains elusive. 
Generally speaking, there are existing at least two difficulties in resolving this issue. First, until now, only two-particle 
entanglement is well explored. How to quantify the multi-particle entanglement is not clear. Second, QPT closely 
relates to the notorious many-body problems, which is almost intractable analytically except in some special toy 
models, such as the Ising model and the XY model. Until now, the only effective and accurate way to deal with QPT 
in critical region is the density- matrix RG method |32j |. Unfortunately, it is only efficient for one-dimensional cases 
because of the much more complicated boundary conditions for two-dimensional situation. 

In order to investigate the entanglement behavior in QPT for two-dimensional cases, we used the original real-space 
RG technique |3lj . This method contains uncontrollable approximations. But by using bigger block size while at the 
same time comparing some of the obtained results with those already known from other methods to check the errors, 
we can still get reliable conclusions. For the entanglement measure, we used the Zanardie's measure. It is in essence 
a bipartite measure, but since most of the present quantum computation applications are concerning two-particle 
entanglement, it is still very meaningful and informative to carry out the study with it. 

IV. CALCULATION METHOD 

We have shown previously that for a half-filled triangular quantum lattice, there exists a metal-insulator phase 
transition with the ratio of electron repulsion and hopping term to be the tuning parameter |33J, l34l l35| . In the 
following, we will use the same model and technique to explore the scaling properties of the entanglements. 

The model we use is the Hubbard model with the Hamiltonian, 

H = -t J2 [t&Cj. + H.c] 

<i,j>,a 

i i 

where t is the nearest-neighbor hopping term, U is the local repulsive interaction and K — — U /A and l. L is the unit 
operator. c^.(ci CT ) creates (annihilates) an electron with spin a in a Wannier orbital located at site i; the corresponding 
number operator is n icr — cf^c^ and < > denotes the nearest-neighbor pairs. H.c. denotes the Hermitian conjugate. 
The order parameter for metal-insulator transition(MIT) is the charge gap defined by A g = E(N e — 1) + E(N e + 
1) — 2E(N e ), where E(N e ) denotes the lowest energy for a N e — electron system. In our case, N e is equal to the site 
number N s of the lattice. We already know that at the critical point (U/t) c — 12.5, there will be a quantum phase 
transition, i.e. MIT, signature by a jump of A g from Zero to non-zero. Unlike the charge gap calculated from the 
energy levels, the Zanardi measure of the entanglement is defined upon the wave function corresponding to E(N e ) 



3 



instead, as show in Eq. ©. Before showing how to get the entanglement from RG method, let us first briefly review 
the multi-staged real-space RG method. 

The essence of real-space RG method is to map the original Hamiltonian to a new Hamiltonian with much fewer 
freedoms which keeps the physical quantities we are interested in unchanged |3l|. The map can be iterated until the 
final Hamiltonian can be easily handled. The crucial step in this method is how to related the parameters between 
the old and the new Hamiltonians. This can be realized by dividing the original lattice into blocks and then build the 
new Hamiltonian upon blocks, namely regard each block to be an effective site. 

Fig.l shows schematically the hexagonal block structure we used in our calculations and the coupling between 
blocks. For each block, we solve it numerically in the subspaces of 6 electrons with 3 spin up and 3 spin down, 7 
electrons with 4 spin up and 3 spin down, 7 electrons with 3 spin up and 4 spin down, and 8 electrons with 4 spin 
up and 4 spin down. In all the subspaces, we keep the lowest-energy non-degenerate state, which is also required to 
belong to the same irreducible representations of C§ v symmetry group. It should be mentioned here that if the there 
is degeneracy then one possible solution is to average over the degenerated states. The kept states will then be taken 
as the 4 states for an effective site, |0) , |f ) , J.) , |f J.) . If denoting the energies corresponding to the first two states 
by Ei,E2, after some intensive calculations, we can obtain the new Hamiltonian on the effective lattice, which has 
the same structure as the original Hamiltonian, i.e., 

<i,j>,a 

+u'Y J (\-<d(\-<i)+ K 'Y. 1 ^ ( 5 ) 

i i 

where the prime ' denotes the operator action upon the block states and 

t' = 3A 2 t, (6) 
U' = 2{E 1 -E 3 ), (7) 
K' = {E x +E 3 )/2, (8) 

with A =< — cr, a\' p c^ {p)ry \ p — a >' p =< <j\' p c^ p) JO >' p , where a denotes f or J.. The above equations are the so-called RG 
flow equations. Usually, they are iterated until we get the fixed point. But besides the fixed points, we can get more 
if we keep track of the results from each RG iteration. For example, if we stop the RG flow at the first iteration, then 
obtained t' and U' will be parameters of an effective 7-site hexagonal block mapped from the original system of 7 2 sites. 
Because the hexagonal block Hamiltonian can be solved numerically, then the physical quantities we have interest for 

a system containing 49 sites can be obtained. In this way, we can study the system of the size 7 2 , 7 3 , 7 4 , 7 5 , 7 6 

This is the so-called multi-stage real-space RG method, which is well adapted to start the finite-size scaling analysis. 
For entanglement, because at the ith iteration, the effective site involves 7 l starting sites, what we are calculating by 
this procedure is equivalent to get the entanglement between the L-site block and the surrounding seven L-site blocks. 
As L increases with the RG flow, the scaling of the entanglement with L comes out readily. 



V. RESULTS AND DISCUSSIONS 



In Fig. 2 we show the results of the entanglement as a function of (U/t) for different system sizes. The crossing 
point in Fig. 2(a) represents the critical value of U/t, which is shown to be (U/t) c — 12.5. It is exactly equal to the 
critical value for metal-insulator transition (MIT) when the same order parameter U/t is used |33l l34l l35j | . In Fig. 2(b), 
with proper scaling, all the curves in Fig. 2(a) collapse onto one curve, which can be expressed as 

E 7 = f{qN^) (9) 

where q = U/t — (U/t) c measures the deviation distance of the system away from the critical state. By the one- 
parameter scaling theory, near the phase transition point, E 7 can written as 

E 7 = q VE f{j) (10) 

where £ = q~ v is the correlation length of the system w here the critical exponent v. Hence, 



E 7 = qy-f(Lqn=q VE f(N^q), 



(11) 
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where we used N — I? for the two-dimensional systems. Now we can have, 

VE = 0,i/= 1 (12) 

It is interesting to note that here we have acquired the same v as in the study of MIT. This shows the great 
consistency of the results since the critical exponent v is only dependent upon the inherent symmetry and dimension 
of the investigated system. 

Another significance of the results lies in the finding that in the metal state, the system is highly entangled with 
Ej = 2 while in the insulating state, the system is partly entangled with Ej — 1. Because the reduced density matrix 
pi is four dimensional, the maximally entangled state can be written as 



P7 



(13) 



(14) 
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in the basis |0) , ||) , \l) , |f j). The related entanglement is Ej — — \ l°g2 \ — 2- Unlike the metal state, in which 

»=i 

the sites can have equal probability in any of the four basis states, the insulating states should be expected to have 
electrons showing less mobility. From the calculations, we know that in the insulating state, 

P7 = 2 1 

2 

V °/ 

which means that the central site has equal probability to be in and no occupation in |0), |TJ.). This well 

3 

coincides with the above expectations. The corresponding entanglement is then E 7 = — ^°§>2 \ = 1- 

All the above discussions are confined to the entanglement between the central site and its surrounding sites. 
Because the central site is a very special site showing the highest symmetry in the block, one may wonder what can 
happen to the neighbor sites, for example, the entanglement between site 1 and the rest 6 sites. To answer this 
question, the same calculations are conducted and the results are presented in Fig. 3. Following the same procedures 
as for site 7, we can get the transition point and the critical exponents. It is no surprise that the same values are 
obtained. The only difference is that in the metal state, the maximal entanglement is a little less than 2 and the 
minimal one is a little less than 1. This can be explained by the asymmetric position of site 1 in the block. 

It should be mentioned that the calculated entanglement here has a corresponding critical exponent i/e = 0. This 
means that the entanglement is constant at the critical point over all sizes of the system, which can be well seen from 
Fig. 2(a) and Fig. 3(a) also. But it is not a constant over all values of U/t. There is an abrupt jump across the critical 
point as L — > oo, which can be easily seen in Fig. 2(b) and Fig.(3)b. If we divide the regime of the order parameter 
into non-critical regime and critical regime, the results can be summarized as follows. 1) In the non-critical regime, 
i.e. U/t is away from (U/t) c , as L increases, the entanglement will saturate finally onto two different values depending 
upon the sign of U/t — (U/t) c . 2) At the critical point, the entanglement is actually a constant independent of the size 
L. These properties are qualitatively different from the single-site entanglement discussed by Osborne [28|, where the 
entanglement with Zanardi's measure increases from zero to the maximum at the critical point and then decreases 
again to zero as the order parameter 7 for XY mode is tuned. The peculiar properties of the entanglement we have 
found here can be of potential interest to make an effective ideal entanglement "switch". For example, with seven 
blocks of quantum dots on triangular lattice, the entanglement among the blocks can be regulated as " 0" or " 1" 
almost immediately once the tuning parameter y crosses the critical point 12.5. The switch errors will depend upon 
the size of the blocks. Since it has already been a well-developed technique to change y f° r quantum dot lattice, 
the above scheme should be well workable. To remove the special confinement we have made upon the calculated 
entanglement, namely only the entanglement of site 1 and site 7 with the rest sites are considered, in the following, 
we will prove that what we are switching here is actually the average pairwise entanglement between the seven blocks. 

Let Eij denotes the entanglement between ith and jth site. Then the total two-site entanglement E to t can be 
obtained, i.e. 

777 
Etot — ^ En + ^ E 2 i + ^ E 3i 

i=2 i=3 i=4 

7 7 

+ Y / E il + Y, E ^ + E e7- (15) 



i—5 i—6 
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From the symmetry of the sites, we know that 



E13 = E24 = -E35 = E46 , En — E25 = E36 
En — E 2 7 — E 37 — E 47 = E 57 = E 67 . 



(16) 



By substituting the above equalities into Eq. (|15|l . we have, 



6E 17 + 3{2E 12 + 2E 13 + E 14 ). 



(17) 



Because 



E 7 — 7E\ 7 , Ei — 2E\2 + 2i?i3 + E 4 



(18) 



We finally have E t 



'tot — 



6E 7 + 3Ei . The average 2-site entanglement is E t 



average — 



^f- = \{2E 7 + E 1 ). Because 



E 7 = f(qN-2), 
E 1 =g{qNi), 



(19) 



then we should have E average — h(qN ' 5 ). This tells us that the average pairwise entanglement also has the properties 
shown in Fig. 2 and Fig. 3. 

Until now, almost all the finite-scaling work about the entanglements are done in one-dimensional cases, such as 
Ising model with transverse magnetic field or XY model. And different scaling properties of the entanglement are 
found there. In |3(|. the one-dimensional Harper Hamiltonian is investigated. The average concurrence shows dif- 
ferent scaling properties in the metal and insulating states. In |29j. Vidal calculated the entanglement with Zanardi 
measure between a L-site block with all the rest spin chains using XY model. And at the critical point, the entan- 
glement demonstrates a logarithm dependence over the size L. Similar results are also found by Osterloh 30] in the 
concurrence derivative over the order parameter A. None of the above work has reported the switching properties of 
the entanglement scaling over size. Hence it seems that the dimension of the investigated systems seem to play an 
important role also in the entanglement scaling, just as expected from QPT viewpoint. 



In summary, by using a multi-stage real-space renormalization group method, the scaling properties of the entan- 
glement with the Zanardi's measure over the block size when QPT happens has been investigated in details. The 
critical exponent v — 1 has been found, which coincides well with our previous work in studying a quite different 
physical property, the charge gap. When the block size L — > 00, the entanglement shows an abrupt change when the 
tuning parameter crosses the phase transition point. This property might be well applied to make an "entanglement 
switch" . 



VI. SUMMARY 



Acknowledgments 



We would like to acknowledge the financial support of the Office of Naval Research(N00014-97-2-0192) 



6 



P. Benioff, J. Stat. Phys. 563, 22 (1980); Phys. Rev. Lett. 48, 1581 (1982). 
R. P. Feynman, Int. J. Theor. Phys. 21, 467 (1982). 

D. Deutsch, Proc. R. Soc. A 400, 97 (1985); Proc. R. Soc. A 425, 73 (1989). 
C. H. Bennett and D. P. DiVincenzo, Nature 404, 247 (2000). 

C. Monroe, Nature 416, 238 (2002). 

L. M. K. Vandersypen, M. Steffen, G. Breyta, C. S. Yannoni, M. H. Sherwood and I. L. Chuang, Nature 414, 883 (2001). 

D. Bouwmeester, J. Pan, K. Mattle, M. Eibl, H. Weinfurter and A. Zeilinger, Nature 390,575 (1997). 
T. Jennewein, C. Simon, G. Weihs, H. Weinfurter and A. Zeilinger, Phys. Rev. Lett. 84, 4729 (2000). 
K. Mattle, H. Weinfurter, P. G.Kwiat and A. Zeilinger, Phys. Rev. Lett. 76,4546 (1996). 

M. Lewenstein, D. Bruss; J. I. Cirac, B. Kraus, M. Ku, J. Samsonowicz, A. Sanpera and R. Tarrach, J. Mod. Opt. 47, 
2841 (2000). 

M. Horodecki, P. Horodecki, R. Horodecki, e-print, quant-ph/0109124 

M. J. Donald, Michal Horodecki, Oliver Rudolph, e-print, quant-ph/0105017 



B. M. Terhal, e-print, quant-ph/0101032 
D. Bruss, J. Math. Phys. 43, 4237 (2002). 

M. A. Nielsen, C. M. Dawson, J. L. Dodd, A. Gilchrist, D. Mortimer, T. J. Osborne, M. J. Bremner, A. W. Harrow and 
A. Hines, Phys. Rev. A 67, 052301 (2003). 

A. Lamas-Linares, J. C. Howell and D. Bouwmeester, Nature 412, 887 (2001). 

M. Eibl, S. Gaertner, M. Bourennane, C. Kurtsiefer, M. Zukowski and H. Weinfurter, Phys. Rev. Lett. 90, 200403 (2003). 
A. J. Berkley, H. Xu, R. C. Ramos, M. A. Gubrud, F. W. Strauch, P. R. Johnson, J. R. Anderson, A. J. Dragt, C. J. Lobb, 

F. C. Wellstood, Nature 300, 1548 (2003). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

J. Schliemann, J. Ignacio Cirac, M. Kus, M. Lewenstei n and D. Loss, Phys . Rev. A64, 022303 (2001). 

J. Schliemann, D. Loss and A. H. MacDonald, e-print, cond-mat/0009083 

Y. S. Li, B. Zeng, X. S. Liu and G. L. Long, Phys. Rev. A 64, 054302 (2001). 

R. Paskauskas and L. You, Phys. Rev. A 64, 042310 (2001). 

J. R. Gittings and A. J. Fisher, Phys. Rev. A 66, 032305 (2002). 

P. Zanardi, Phys. Rev. A 65, 042101 (2002). 

S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 1999). 

T. J. Osborne and M. A. Nielsen, e-print, |quant -ph/01090'24| 

T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 03 2110 (2002). 

G. vidal, J. I. Latorre, E. Rico and A. Kitaev, e-print, quant-ph/02 11074 
A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 416, 608 (2002). 

Real-Space Renormalization (Topics in Current Physics), Editors: T. W. Burkhardt and J. M. J. van Leeuwen, Springer- 
Verlag, 1982. 

S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

J. X. Wang, S. Kais and R. D. Levine, Int. J. Mol. Sci. 3, 4 (2002). 
J. X. Wang, S. Kais, Phys. Rev. B 66, 081101(R) (2002). 

J. X. Wang, S. Kais, F. Remade and R. D. Levine, J. Chem. Phys. B 106, 12847(2002). 
A. Lakshminarayan and V. Subrahmanyam, Phys. Rev. A 67, 052304 (2003). 



7 



1 2 




7 2 sites 



7 n sites 



FIG. 1: Schematic diagram of the triangular lattice with hexagonal blocks. The hexagonal block is used for the renormalization 
process. The iteration direction is shown by the arrow. The dotted lines represent the interblock interactions and solid line 
intrablock ones. 
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FIG. 2: Variations of the entanglement between site 7 and the neighboring sites (See Fig.l) against the ratio of on-site electron 
interaction U to the hopping term t for different system size, i.e. the number of sites: 7(solid circle), 7 2 (open circle), 7 3 (solid 
square), 7 4 (open square) 7 5 (solid down triangle), 7 6 (open down triangle), 7 7 (solid up triangle), 7 8 (open up triangle), 7 9 (solid 
diamond) 
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FIG. 3: Variations of the entanglement between site 7 and the neighboring sites (See Fig.l) against the ratio of on-site electron 
interaction U to the hopping term t for different system size, i.e. the number of sites: 7(solid circle), 7 2 (open circle), 7 3 (solid 
square), 7 4 (open square) 7 5 (solid down triangle), 7 6 (open down triangle), 7 7 (solid up triangle), 7 8 (open up triangle), 7 9 (solid 
diamond) 




FIG. 4: The same as in Fig.2, but for the entanglement between site 1 and the rest 6 sites of the block. 
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FIG. 5: The same as in Fig.2, but for the entanglement between site 1 and the rest 6 sites of the block. 



